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Abstract 

We present a factorization formula for valence quark distributions in a hadron in x ^ 1 limit. 
For the example of pion, we arrive at the form of factorization by analyzing momentum flow in the 
leading and high-order Feynman diagrams. The result confirms the well-known 1 — x scaling rule to 
all orders in perturbation theory, providing the non-perturbative matrix elements for the infrared- 
divergence factors. We comment on re-summation of perturbative single and double logarithms in 
1 — 
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Three decades ago, Brodsky, Farrar and others made a perturbative quantum chro- 
modyanmics (pQCD) prediction about the power behavior of parton distributions as x — 1, 
i.e., the density of quarks carryingalmost all the longitudinal (plus) momentum of a hadron 
participating in hard scattering [H 0, S 0] • Through calculating the leading pQCD dia- 
grams, they showed that the valence quark distribution goes like, for example, (1 — x)^ in 
the pion and (1 — x)^ in the nucleon. The basic argument is that when the valence quark 
carries nearly all the plus momentum, the relevant QCD configurations in hadron wave 
functions must be far off-shell and hence are amenable to pQCD treatment. The result has 
since been generalized to the sea quarks, gluons, helicity-dependent distributions 0, and 
lately to generalized parton distributions p. The above power law is consistent with the 
large-momentum behavior of the elastic form factors through the so-called Drell-Yan-West 
relation ^, ^, ^J. 

Despite its elegance and usefulness, a rigorous derivation of the scaling rule in pQCD 
has not, to the authors' knowledge, been established in the literature. Indeed, if one follows 
the original calculations, the scaling term has a power-like infrared-divergent coefficient 
J d^kx/ {k\)'^ for the pion, and similarly for the nucleon. One might argue that this di- 
vergence will be tamed by non-perturbative QCD effects and hence does not affect power 
counting. A more satisfactory solution, however, is to formulate a factorization theorem 
explicitly separating the long distance contributions from short distance ones. QCD factor- 
ization will establish the scaling rule through a systematic QCD power counting, and hence 
validate the result to all orders in perturbation theory. It will also provide well-defined 
non-perturbative matrix elements to absorb infrared divergences. Finally, factorization will 
provide a convenient tool to re-sum possible large Sudakov double logarithms present in 
X — 1 limit. 

The goal of this paper is to establish a factorization formula for the valence quark distri- 
butions in a hadron in the large-x limit. We consider the example of a pion, although similar 
studies apply to other hadrons such as the nucleon. The factorization arises from a detailed 
analysis of momentum fiow in the lowest-order pQCD diagrams and their generalization to 
all orders. Without an explicit calculation, an infrared power counting allows a straightfor- 
ward determination of the 1 — x power behavior. The factorized expression contains hard 
parts, initial and final state jets, and a soft function, as easily seen from the leading 1 — x 
structure of reduced diagrams. We introduce a non-perturbative matrix element for the soft 
factor which absorbs the perturbative infrared divergence mentioned earlier. 

The large-x behavior of the valence quark distributions in a pion was obtained from 
calculating the lowest-order Feynman diagrams, one of which is reproduced in Fig. la. 
Let us choose the light-cone system of coordinates defined by two light-cone vectors = 
A(l, 0,0,1) and nf" = (1, 0, 0, -1)/(2A) {p"^ = n"^ = 0, p ■ n = 1). Any vector k^" can be 
expanded in terms of and n'^, as k'^ = k-np^ + k-pn^ + k^. We choose the pion momentum 
mainly along the p-direction (p"*" taken as large), = p^ + m^/2n^, ignoring the m^-term 
except when it is needed to regularize infrared divergences. As labelled in the figure, the 
incoming quark carries momentum Xip^ and anti-quark momentum (1 — Xi)p^ = Xip^. The 
anti-quark going through the final state cut (shown by the dashed line) has momentum k^ 
integrated over. 

Instead of doing an explicit calculation to reproduce the familiar result g(x) ~ /^(l — 
x)^ J (Pkx/lk'^)'^, let us analyze the momentum fiow of the loop integral over k^. In the 
X — ^ 1 limit, the p (or plus) component of k is constrained to (1 — x)p and is soft, ~ A. If 
the transverse momentum is on the order of (1 — xY^'^p^, the n component of k must 
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FIG. 1: a). A lowest-order Feynman diagram contributing to the leading-1 — x power of valence 
quark distribution in a pion. b). A reduced diagram in which the final state anti-quark is a jet 
along the n-direction. c). An alternative reduced diagram in which the final state anti-quark is 
soft. 



be very large, k~ = k']_/2{l — x)p^ ~ p^. Then the anti-quark hne is jet-hke, and going 
along the n-direction. The corresponding reduced diagram is shown in Fig. lb. If, on the 
other hand, k± is soft (1 —x)p'^, k~ also remains soft. The anti-quark line is now a soft line. 
The corresponding reduced diagram is shown in Fig. Ic. Using an infrared power counting 
generalizable to all orders in perturbation theory, one can obtain g(a;) ~ (1 — xY . We refer 
the reader to well-known references for the details of infrared power counting 



If the anti-quark is coUinear to the ra-direction (Fig. lb), the infrared power of the diagram 
is 0, obtained from the final-state quark line as a jet, J d^k6{k'^)6{k^ — (1 — x)p'^). The 
polarization sum of the final-state quark is proportional to k~ ~ p"*" and hence does not 
contribute any soft power. However, an explicit calculation also finds a factor k'j_ in the hard 
part. This suppression factor is not included in ordinary power counting, and is resulting 
from the spin structure of the diagram. Taking into account this additional contribution, the 
actual infrared scaling of Fig. lb is A. Infrared factors of a Feynman diagram come from both 
1 — x and transverse momentum. The dimensionless quark distribution is proportional to 
from the pion wave function. Thus the transverse-momentum integral must be proportional 
to / dPk^/ {k\Y . Hence, the quark distribution has (1 — a;)-related infrared power 2, or 
is proportional to (1 — x)^. The same counting applies to all other leading-order diagrams 
which are not shown explicitly. 

When the final-state quark is soft, the gluon and internal quark lines become coUinear 
to the pion momentum (Fig.lc). The denominators in the pion jets contribute A~'^. Two 
collinear-quark-gluon couplings contribute soft-power A^: Usually a coupling contributed 
by a coUinear transverse-momentum is counted as A^/^. In the present case, however, 
the transverse momentum in the pion jets originates from the soft quark and thus must be 
treated as the soft scale A. The final-state soft quark contributes A^. So the total soft power 
is 0. After factoring the transverse-momentum integral, one has a (1 — a;)-related soft-power 
2, or the pion distribution goes like (1 — x)^, just like in the first region. Fig. lb. 

The above consideration and result are generic. For a general Feynman diagram of an 
arbitrary order contributing to parton distributions at x — >■ 1, the momentum integrations 
contain many different regions (pinched surfaces) 0, 11 1, each of which can be represented 
by a product of the various parts: the two incoming jets along the p-direction {Jl,r), the 
out-going jets along the n-direction (Ja), and two hard parts {Hl^h), each on the different 
sides of the cuts connecting the p and n jets, and finally a soft part connecting aU jets and 
hard parts, as shown in Fig. 2a. Similar factorization for the quark distribution in a quark 
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FIG. 2: a). A general reduced diagram for parton distributions at x — > 1. b). A leading region 
where the anti-quark goes through the soft part and the final state jet containing the eikonal line 
and longitudinally polarized gluons. b). Another leading region where the anti-quark is part of 
the final-state jet. There is no soft factor here because the initial and final state jets are all color 
neutral. 



has been considered in l^, 



Let us use pj to denote the number of coUinear quarks or gluons with-physical-polarization 
hues entering the hard part from jet J; use Ij to represent the number of colhnear gluons 
of longitudinal polarization through a similar attachment; use Ej'^ to denote the number 
of soft bosons or fermions connecting to the left or right hard part; use -E'^'^ ^ to label the 
number of soft bosons of fermions connecting to the left or right hard part; and finally use 
Vj to label the number of three-point vertices in the jet, and sj the number of soft gluons 
with scalar polarization and soft quarks attaching to the jet. Then the soft part has infrared 
power, 

= EjL + Ejji + E\p + E]jj^ + E]jj^ 

+ \[E\^ + E\^ + E\p + E],^ + E],^) . (1) 

On the other hand, the infrared power associated with each jet is 

= \{PJ -sj-E)-Efj-Ij) + ]^ [sj + lj- e {sj + lj- v^'^) . (2) 
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Combining the above results, one finds 

u^{G) > ^ {pj^ + pj^ + pj^) - 2 - + JjJ 

+1 [e]^ + + EX) - \ [sj, + sj, + 

where we have taken into account the constraint that the plus momentum going through 
the final state cut is (1 — x)p ~ A, and the eikonal line is accounted for by pj^ = 2. The 
above result generalizes that for a single quark state jl^ . 

For a pion, /j^ = Ij^ = 2, and the lowest infrared power is —1, obtained when pj^ ^ = 1, 
PJ/3 = 2, E^j^ = Ejj^ = 1, = 0, sj^ = E^, and s = E^^^ + 1. The corresponding 
reduced diagram is shown in Fig. 2b, where the final-state jet contains just the eikonal line 
and longitudinally-polarized gluons, and the soft part contains an anti-quark line. The jets 
going into the hard parts contain a single quark and an arbitrary number of longitudinally 
polarized gluons. We argue, however, that the actual infrared power is when taking into 
account an additional numerator suppression, just as in the leading-order diagram in Fig. 1. 
Observe that the infrared power of an initial-state jet is always an half-integer, which means 
the jet is proportional to at least one transverse momentum, even after all integrals in the 
jet are carried out. However, the only relevant transverse momenta available are those of 
the soft lines. Therefore, this transverse momentum must be counted as A, not A^/^. With 
this extra suppression, the actual infrared power of the reduced diagram in Fig. 2b is 0. 
Factorizing out the transverse-momentum integral, one gets q{x) ~ (1 — xY to all orders in 
perturbation theory. 

The next lowest infrared power, 0, is obtained when pj^ = pj^ = 2, pj^ = 4, all E-f = 
and E'^ = 0, and Ej = sj. The corresponding reduced diagram is shown in Fig. 2c, where 
the final-state jet has an anti-quark quantum number. The soft-gluon radiation is not drawn 
because the initial and final state jets are color-neutral; and thus when summing over all 
final states, the soft factor reduces to 1, i.e., the soft radiation must cancel. Because of the 
angular momentum conservation, there is a numerator suppression factor proportional to 
the transverse-momentum k\ of the ant i- quark jet, as seen in Fig. la. The argument goes as 
follows. Consider deep-inelastic scattering on a pion target. If the final state consists of two 
collinear quark-antiquark jets plus an arbitrary number of longitudinally polarized gluons 
going exactly along the n direction, the scattering amplitude vanishes identically. This is be- 
cause when the initial and final states are collinear and with zero helicity (fermion helicity is 
conserved in gauge theory), scattering with a physical photon of helicity ±1 is forbidden by 
angular momentum conservation. Therefore, there must be a relative transverse-momentum 
between the quark and anti-quark jets, k"^ ~ (1 — x)Q'^. The scattering amplitude is propor- 
tional to k±, and the parton distribution is proportional to k']_. The actual infrared factor 
is then A. Accounting for the transverse-momentum integral, one again gets a scaling law 
(1 — x)^ to all orders in perturbation theory. 
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One can write down a factorization formula for the two reduced diagrams shown in 
diagrams 2b and 2c, after factorizing the colhnear longitudinally polarized gluons from the 
hard parts and the soft gluons from the initial state jets. For example, let us consider the 
case when the anti-quark is soft (Fig. 2b). Going through steps similar to those outlined in 
Tlilll], one finds 



q^{x^l) = HL{p,fx)HR{p,n)JpAp,jj)JpAp,n) 



xE n 



(27r)" (27r)" (27r)" 
X [dy [ dzS^'''^^yp,^^;{ql-^J;{ql:-J;{qn)Jf'\zp,^^;n;^ 



\ rriL niR I 



X^" y.€mr. +yyR,m. + V 5(1 - X - |/ - .) , (4) 



where Hl r represent the hard contributions with collinear longitudinally-polarized gluons 
factorized, Jl^r represent the jet contributions with soft gluons and quarks factorized. The 
factorized soft gluons are now represented by the two eikonal lines £l,r- The soft loop 
momenta to be integrated over connect the soft part and the eikonal lines, the soft part and 
final-state jets. The summations are over all possible cuts through final-state jet Cp and 
soft part Cs- The longitudinal momenta going through the final state cuts are constrained 
by 1 — X. The soft function and jets can be represented by non-perturbative QCD matrix 
elements. 

One can define the gauge-invariant quark-eikonal cross section, 

a^^'^((l - = ^ I ^e^(i-)p+J/- 5^(0|^„(r)|^)7+(^l^a(0)|0) , (5) 

n,a 

where \I/(i/~) is defined as 

^(y-) = / dy+U„ioo;y-)UpiO+,y-;y+,y-)i^iy+,y-) , (6) 



with Un{oo;y ) = Pexp (^—ig dXn ■ A{\n)^, etc., and P in Up is the pion momentum. 

0"''''^ has a similar factorization as the parton distribution, except the pion jets are replaced 
by eikonal jets. Therefore, a; ^ 1 quark distribution can be expressed as, 

q{x) = Hl{p, Ij)Hr{p, fi)JL{p, I^)Jr{p, - x)p, /i) . (7) 

The entire x-dependence can now be found from the quark-eikonal cross section. This result 
is ver y si milar to the factorization of the quark distribution in a single quark in the same 
limit llJ, ll^. It is easy to check that at tree level, one reproduces the result from Fig. 
1. 

Apart from the simple power behavior considered so far, perturbative calculations also 
generate single and double logarithmic dependence in 1 — x. The single logarithmic depen- 
dence can be summed to yield a power correction (1 — x)'^^^\ where v4(/i) is renormalization- 
scale fj, dependent jl^. The evolution in can be calculated in perturbation theory, and 
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the absolute magnitude of can only be determined through non-perturbative methods. 
Soft-gluon radiations from charged lines pro duce large double logarithmic corrections of type 
a" ln'"(l — x) with m < 2n — 1 [l^, IfiljlTl Esl . K^ . These corrections depend on a precise 
definition of the parton distribution |l2j| and infrared regularizations. These issues will be 
discussed in a forth-coming publication. 

To summarize, we have presented a factorization formula for the valence quark distri- 
bution in the pion valid to the leading power in (1 — x) and to all orders in perturbation 
theory. Through infrared power counting, one obtains the power dependence in (1 — x). 
The factorization theorem, however, allows one to calculate the coefficients of the power law 
non-perturbatively. It also allows one to re-sum large Sudakov double logarithms associated 
with soft radiations. 
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